Thermal and Dynamical Properties of the Two-band Hubbard Model 

Compared with FeSi 
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We study the two-band Hubbard model introduced by Fu and Doniach as a model for FeSi which 
is suggested to be a Kondo insulator. Using the self-consistent second-order perturbation theory 
combined with the local approximation which becomes exact in the limit of infinite dimensions, we 
calculate the specific heat, the spin susceptibility and the dynamical conductivity and point out 
that the reduction of the energy gap due to correlation is not significant in contrast to the previous 
calculation. It is also demonstrated that the gap at low temperatures in the optical conductivity is 
filled up at a rather low temperature than the gap size, which is consistent with the experiment. 



I. INTRODUCTION 

For long time, FeSi has attracted interests of many re- 
searchers. At low temperatures, FeSi is a nonmagnetic 
semiconductor with a narrow band gap of about 700K, 
and becomes magnetic and metallic when the tempera- 
ture rises. Jaccarino, et al. JO reported an anomalous 
decrease of magnetic susceptibility at low temperatures, 
which could not be explained from conventional band 
theories, unless the vanishingly small band width was as- 
sumed. Moreover, it was found by experiments that there 
was no magnetic ordering at low temperatures. Q Taka- 
hashi and Moriya H applied the spin-fluctuation theory 
to the nearly ferromagnetic semiconductor model, and 
explained the temperature variation of the specific heat 
and the susceptibility. Their idea of thermally induced 
local moment was confirmed later by the neutron scat- 
tering experiment. Q 

On the other hand, several materials with strongly 
correlated f-band were found to exhibit an energy gap 
and become insulators at low temperatures. They are 
called Kondo insulators and have been studied inten- 
sively. Typical examples are YbBi 2 , Cc3Bi 4 Pt3 || 
and many others. Aeppli and Fisk Hf suggested that 
FeSi may be viewed as a Kondo insulator in d-electron 
systems, a band insulator with moderate electron corre- 
lation. Schlesinger, et al. || measured the optical con- 
ductivity and found that the gap (~ 700K) at low tem- 
peratures is filled up at a temperature (~ 250K) lower 
than the gap size. Therefore, they concluded that FeSi 
could not be understood with simple band models with- 
out correlation effect. Low temperature properties were 
elucidated further, ||,[l0| and the importance of the cor- 
relation effects was underlined. jll| 

The band calculation by Mattheis and Hamann |l^] 
gave the indirect gap of about O.lleV, which is slightly 
larger than the observed one. Fu, et al. |L3| also calcu- 
lated the indirect gap of about O.leV. The band width, 
however, was much larger than that assumed in the anal- 
ysis of ref . . 

Recent ab initio band calculation with spin-orbit inter- 



action, Jl4|,[l5| however, has given a rather small gap size 
(indirect gap - 0.02-0.03eV, direct gap ~ 0.046-0.08eV). 
It is well known that the local density approximation un- 
derestimates the gap of insulators, (l6) and the precise 
estimate of the small gap size is a difficult task. There- 
fore, it seems reasonable to consider that the gap size 
obtained by the band calculations are of the order of the 
observed value and no large reduction of the gap due to 
correlation seems necessary. 

Some of the low temperature properties were calcu- 
lated from the band calculations |l7j , and the agreement 
with experiments were obtained. Temperature depen- 
dence of the optical conductivity was calculated from the 
band calculations also, [ ^3|Jl^1 but the experimentally ob- 
served behavior || was not reproduced. 

On the other hand, high-resolution photoemission 
measurement was done by Saitoh et al. They esti- 
mated the self-energy of the d-electrons from the differ- 
ence between the photoemission spectra and the previ- 
ous band calculation, and demonstrated a possible 
existence of the strong renormalization of the electronic 
states at low energies. 

In the above band calculations, minimum gap is given 
by the indirect gap, but it is close to a direct gap, and 
the states near the Fermi level mainly consist of the d 
electrons. Based on these, Fu and Doniach proposed 
a two-band Hubbard model. In a non-interacting case, 
it has a simple direct gap originating from a hybridiza- 
tion between the two bands, and the Coulomb repulsion 
was assumed to act only within each band. The inter- 
band Coulomb repulsion and the exchange interaction 
were simply neglected. Therefore, this is the simplest 
model for FeSi, but may keep essential features near the 
Fermi level. They studied many body effects using the 
perturbation theory with respect to the Coulomb inter- 
action between d-electrons. But their calculations seem 
to include some errors in the treatment of the proper 
symmetry of the self-energies. 

The purpose of the present paper is to investigate their 
model correctly and to study the effects of correlation on 
the electronic structures and the physical properties. The 
correlation effects are not included in the spin fluctuation 
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theory @] at low temperatures. We calculate the temper- 
ature dependence of the density of states, magnetic sus- 
ceptibility, the specific heat and the dynamical conduc- 
tivity using the self-consistent second order perturbation 
theory (SCSOPT) Q together with the local approxi- 
mation, which becomes exact in the limit of large spatial 
dimensions, plj] Our basic formulation is the same as 
the previous one fll9|| , but the obtained results are much 
different. Our main conclusions are: (1) the gap is not 
much renormalized by the correlation, and (2) the strong 
temperature dependence of the density of states due to 
correlation results in the rapid collapse of the gap in the 
optical conductivity. Although the model used here is a 
simplified one, these results can explain the experiments 
semi-quantitatively. Similar calculations were performed 
on the periodic Anderson model, p2] , p3t and the rapid 
gap-filling was found. Comparison of the present results 
with these calculations are also of interest. 



U(n a i + n a i)/2 in Ti a cancels with the Ed term. In 
the case of a finite magnetic field, we must calculate n atT 
self-consistently. 

In SCSOPT, is calculated as 



e2» = w 



de 1 de 2 de 3 pj(e 1 )p^ b (e2)pj{e 3 ) 



x /(-ei)/(e 2 )/(e 3 ) + /(gi)/(-g 2 )/(-£ 3 ) ^ (g) 
uj + E\ — s 2 — £3 + i<5 



where 



and 



Pab{") = ImG£ 6 (w + tf) 

7T 



(4) 
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II. MODEL AND FORMULATION 



Here, N is the number of sites. These equations are 
converted into MM 



The two-band Hubbard model introduced by Fu and 
Doniach jl9| is written as 

ija 



E$» = -it/ 2 / dre- [ S a T(-T)^ 6 (r)A- 6 CT (r) 



(1) 



where cj acr (ci Qcr ) creates (destroys) an electron on site i 
in band a = 1, 2 with spin cr. The opposite signs in the 
hopping terms of the two bands provide a simple direct 
gap when they are hybridized via v. The magnetic field 
is denoted by h. Green's functions in the k-space can be 
represented in the matrix form: 



_ ^w-et-SfiOc.w) 



-«-Ef 2 (k,w) 

-£ k -E£ 2 (k, W ) 



(2) 



We 



where £k denotes the Fourier transformation of tij. 
choose tu = 0. 

To investigate the effects of electron correlation in this 
model, we adopt SCSOPT combined with the local ap- 
proximation, pfj The self-energy consists of the first and 
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the second order terms, E^ fc (w) 

where Ea = U(n a -\ + ?Vl)/ 2 + Ua(n a ^ - n a {)/2, and 
n atT denotes the number of electrons in band a with spin 
a. Wc implicitly introduce d- level energy Ed = —{7/2 
to keep particle-hole symmetry. Then the first term 



where 



+A2(-r)Bl b {r)B-J{r) 

/OO 
dee-™pl b {s)f{e), 
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B° ab {r) = / dee-™ P l b {e)f{-e). 
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The basic formulation mentioned above is the same as in 
rcf. p9| , but further analysis seems different from ours 
indicated below. 

In the particle-hole symmetric case, each element of 
the Green's functions is calculated as follows: 

G?iM = G&M = ^-^W.H), (8) 



GJaH = GJiM = ^^F(C.H), (9) 



E^H^E^H and E? 2 (w) = E^w), (10) 



where 



1 



(11) 



To carry out k-summation, we assume the semi- 
circular shape of the density of states ( D(e) = 
(2/7rW)yl — e 2 /W 2 ) for the unperturbed band without 
hybridization. Then we obtain 
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F{z) = —( z -^/z 2 - W 2 



(12) 



Though Fu and Doniach assumed the Gaussian form to 
the density of states, the difference in the shape of D(e) 
is not essential for the purpose of our calculation. 

In the non-interacting case, the self-energy vanishes, 
so that we obtain the density of states as 

P11M = P22M 



\/uj 2 — v 2 



D{^Jlo 2 -v 2 )9(\uj\ -v), 



P12M = P21M 

usgn(o;) 



y/ui 2 — v 2 



D(Vtu 2 -v 2 )6{\u\ -v). (13) 



Note that Pi 2 (oj) is an odd function. 



III. CALCULATION OF PHYSICAL 
QUANTITIES 

Starting from the equation of motion, p4| we obtain 
the following expression for the total energy per site: 

6 = 9 E / ^/M ["KM + P22M} 

\ E { £ kP?i(k, w) + (-£ k )^ 2 (k, w)}] , (14) 
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where 



^ 6 (k,w) = --ImGE 6 (k,w + W). (15) 
For [7 = 0, the above equation reduces to the trivial form, 

/oo 
dw/Mw^H +/&(«)]. (16) 
-00 



The k-sums in the brackets in eq.(p^) are calculated as 
follows: 



l^ekG-^k^) = I^(- £k )G£ 2 (k,a;) 

= -l + Ca(w)F(C tr (w)). (17) 
The specific heat is given by 

dt 



dT 



(18) 



This differentiation is carried out numerically. 
The magnetic moment per site is given by 



m(T, h) 
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) 

+ pl 2 (u;)-p l 22 (u;)}. (19) 



The magnetic susceptibility is calculated by x{T) — 
ra(T, h)/h at sufficiently small h. Furthermore, using the 
fact that the vertex correction in the dynamical conduc- 
tivity drops out in d — 00 theory, |2l| we also obtain the 
expression for the dynamical conductivity in Fu-Doniach 
model as 

Rea^T) = 2^ f dedvDje) f {v) ~ & + " } 

x \pn (e, ( £ : ^ + w) - P12 (e. f)Pia {e,v + L0% (20) 

where 

pJ 6 (e,w) = -^X;*(ek-e),t&(k,w). (21) 

k 

For [7 = 0, h = this equation becomes 
,.2 



7r/> u (w/2)_5a(l-2/(u;/2)) 



+47ro"(w) / oVD(\/ 



(22) 



The first and second terms describe the interband and 
intraband (Drude) contributions, respectively. 



IV. RESULTS 

We take the unit of energy as W = 1 and choose v = 
0.125, which yields the gap of 2v = 0.25 for U = 0. 

In Figs. 1(a) and 1(b), we show the self-energies. 
Though the diagonal parts have rather similar shapes to 
those of the single-band Hubbard model, the off-diagonal 
parts have the opposite parity compared to the diago- 
nal one. This can be understood from the odd parity 
of (0i2 (w) in eq.fO), since Im£^ b has the same parity as 
Pah(w) (see eq.(g)). In ref. fl9fl , both ImSn and Im£i2 
are treated as even functions, which implies that both 
ReSn and Re£i 2 have odd parity and vanish at e = 
in their calculation. In our proper treatment, ReS 12 is 
an even function and has a positive value at low ener- 
gies and is added to the mixing v in G a b, so that the 
reduction of the gap due to the renormalization factor 
z = (1 — dY* aa (u)/du!\ UJ= o)~ 1 is almost compensated by 
this increase of the effective mixing (see below). 
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FIG. 1. The real and imaginary parts of the self-energies 
Sab(w) are shown for U = 2, v = 0.125 and T = 0. The 
full and the broken lines indicate the diagonal (En) and the 
off-diagonal (E12) parts, respectively. 





FIG. 2. (a) The quasi-particle density of states pn(w) at 
T — is shown for U =0, 1, 2 and 3. The insert indicates 
pi2(w) for T = and U — 2. (b) pn(w) for (7 = 2 at various 
temperatures are shown. 

In Fig. 2(a) and 2(b), the quasi-particle density of 
states obtained with the above self-energies are shown. 
In the present model, the density of states for U = has 
a gap of width 2v that has sharp edges and diverges at 
\oj\ = v. So we take 8 finite (8 = 10 -4 except in <j{uif) in 
lu + iS in order to make numerical calculation easy. It is 
seen in Fig. 2 (a) that the gap size is not significantly re- 
duced when U increases. The reason is mentioned above. 
This result is in strong contrast with the periodic Ander- 
son model, in which the gap is strongly reduced by the 
Coulomb repulsion between f-electrons. |22| , p5| ] 

We have also investigated a variation of the present 
model, where the band energies £k and — £k in the band 
1 and 2 are replaced with and — be^, respectively, and 
the Coulomb repulsion U in band 1 is replaced with bll. 
If we set b = 0, the band 2 becomes dispersionless and the 
band b becomes free, so we obtain the periodic Anderson 
model, in which the gap is known to be reduced due to the 
correlation. We confirmed that the gap is continuously 
reduced when we change from b = 1 to 0. 

The density of states varies with temperature when U 
is finite. In Fig. 2(b), we show the density of states for 
U = 2 at various temperatures. The gap is almost filled 
up at the temperature of the order of the gap. It is filled 
up faster for larger U. 

Figure 3 shows the temperature dependence of the elec- 
tronic specific heat per site. The peak in the U = curve 
splits into lower and higher temperature regions due to 
correlation effect at finite U. This tendency becomes less 
remarkable if we use the density of states at T = to 
calculate Cy(T) at finite temperatures. It indicates the 
importance of the temperature dependence of the density 
of states itself in the calculation of the physical quanti- 
ties. 
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FIG. 3. Temperature dependence of the specific heat is dis- 
played for U —0, 1, 2 and 3. 
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FIG. 4. Temperature dependence of the spin susceptibility 
in unit of (p^is) 2 /W is displayed for U —0, 1, 2 and 3. The 
circles indicate the experimental result [1] (see the text in V). 

The temperature dependence of the spin susceptibility 
per site for different magnitude of U are displayed in 
Fig. 4. As U increases, x(^) is largely enhanced in the 
intermediate temperature region, and the peak positions 
shift to lower temperatures, similar to the specific heat. 

In Fig. 5, the dynamical conductivity for U — 2 at var- 
ious temperatures is shown. At T = 0, <j(uj) has a gap of 
about 0.23, although the gap edge is rounded due to the 
finite S = 2 x 10~ 3 . The sharp peak at the gap edge re- 
flects the peaks in the density of states. a(uj) varies with 
temperature faster than the density of states because it 
consists of the convolution of the two spectral functions 
and the Fermi distribution function. As a result, the gap 
is filled up almost completely at the temperature which 
is about a third of the gap size. If we take a smaller value 
of U, e.g. U = 1, the gap is filled up more slowly, at a 
temperature which is a half of the gap. 
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FIG. 5. The dynamical conductivity for U — 2 is shown 
for T=0, 0.04 and 0.08. The insert shows the comparison 
of the calculation with the experiments at T=20K(circles), 
150K(triangles) and 250K(squares).[8] 

V. DISCUSSIONS 

Although the model used in the present study is over- 
simplified as a model for FeSi, we try to compare the 
obtained results with the experiments. The gap size of 
FeSi is estimated from the activation energy of the mag- 
netic susceptibility and the resistivity. Both yield the gap 
size of about 650-700K. || In the measurement of <j{lo) at 
T =20K, there is a threshold at about 500cm" 1 ~ 700K. 
U Therefore, we first take 700K as the gap at T = 0. 

If we regard the gap size 0.23 for U = 2 in a(uS) to 
be corresponding to the gap of 700K, the peak position 
in x(T') (T ~0.11) corresponds to about 340K, which is 
too low compared the experiment (500K B). If we take 
[7 = 1 with the gap ~ 0.24 in cr(w), the peak position 
in x(T) becomes ~ 0.11 ~ 470K, which is closer to the 
observed one. Recent experiment at a lower temperature 
(T =4K), [^6) however, reported a little larger gap size 
~ 570cm~ 1 ~ 820K. Use of this value results in the peak 
position of x(T) with U = 2 at 400K, which is not much 
different from the experiment (see the circles in Fig. 4). 
The absolute value is fitted by adjusting the value of the 
effective g-factor (g =2.53), which represents the effects 
of the multi-band and the electron number in the real 
system. 

As we have shown in the previous section, the den- 
sity of states in Kondo insulators depend on tempera- 
ture. As a result, the activation energy analysis from 
x{T) or resistivity may underestimate the gap size at 
T = 0. We therefore conclude that U = 2 is more ap- 
propriate to fit both a(u>,T) and x(T) within the present 
model, although refinement of the model and the calcu- 
lation method is definitely needed. 

There is a large ambiguity in the experimental deter- 
mination of the electronic part of the specific heat, 0] 
so that we did not compare the peak position of Cy(T) 
with experiments. 
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In the present paper, we carefully investigated the two- 
band Hubbard model proposed by Fu and Doniach as a 
simplest model for FeSi. The electron correlation was 
taken into account by use of the self-consistent second- 
order perturbation theory. We calculated the temper- 
ature dependence of spin susceptibility, electronic spe- 
cific heat and the dynamical conductivity. The observed 
temperature dependence of x(T) and a(u>,T) were semi- 
quantitatively explained by the present model if we fit the 
optical gap size to the recent experimental one and choose 
U/W = 2. However, <j(lu) in the present calculation has 
a sharp peak at the gap edge, which is not observed in 
experiments. This peak originates from the large peaks 
in the density of states, and would be suppressed if we 
use the density of states obtained from band calculations 
p2|- |l5| ], or if we introduce disorder (^3|. Secondly, the 
filling up of the gap of the measured <j(uj) with increas- 
ing temperature does not seem to be compensated by the 
sufficient decrease in the region right above the gap. A 
possible solution for the case of the periodic Anderson 
model was discussed in plj. Thirdly, the observed er(w) 
has a large tail in the high frequency region, irrespective 
of temperature. To consider these features correctly, one 
needs to start from a more realistic energy band model. 
We like to perform such calculations in our next studies. 
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